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Abstract. 

The problem of chaotic scattering in the presence of direct processes or prompt 
responses is mapped via a transformation to the case of scattering in the absence of 
such processes for non-unitary scattering matrices S. When prompt responses are 
absent, S is uniformly distributed according to its invariant measure in the space of S 
matrices with zero average (S) = 0. When direct processes occur, the distribution of 

5 is non-uniform and is characterised by an average (S) ^ 0. In contrast to the case 
of unitary matrices S, where the invariant measures of S for chaotic scattering with 
and without direct processes are related through the well-known Poisson kernel, we 
show that for non-unitary scattering matrices the invariant measures are related by the 
Poisson kernel squared. Our results are relevant to situations where flux conservation 
is not satisfied, for transport experiments in chaotic systems where gains or losses are 
present, for example in microwave chaotic cavities or graphs, and acoustic or elastic 
resonators. 



PACS numbers: 73.23.-b, 03.65.Nk, 42.25.Bs 
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1. Introduction 

The statistical properties of ensembles of scattering matrices have been studied 
extensively since their introduction in pioneering works in the field of nuclear physics [U 
|2]. More recently, Random Matrix Theory (RMT) techniques have been applied to 
study several statistical properties of electronic transport in mesoscopic systems [HIIUE]- 
Although the spectrum of applications of RMT is quite wide, ranging from atomic nuclei 
to microwave cavities, most of the investigations in RMT have considered systems in 
which flux is conserved. Therefore, the associated n x n scattering matrix S is unitary, 

S& = I n , (1) 

where I n denotes the unit matrix of dimension n, and the dagger means Hermitian 
conjugation. Here n denotes the number of scattering channels. 

When unitarity is the only constraint, the unitary case denoted by (3 = 2 in Dyson's 
scheme applies [6] . If in addition time- reversal symmetry is imposed, then S is also a 
symmetric matrix, 

S = S T , (2) 

a case corresponding to (3 = 1. Here, the superscript T denotes transposition. When 
time-reversal symmetry is present, but spin-rotation symmetry is broken, the S matrix 
is unitary and self-dual, a case denoted by (3 = 4. 
Any scattering matrix S can be decomposed as 

S = UV, (3) 

where U and V are n x n unitary matrices for (3 = 2, while V = U T for (3=1. For 
(3 = 4, V = U R , where U R is the dual matrix of U. 

It has been shown [HI H] that chaotic scattering in the absence of direct processes is 
well described by uniformly distributed S'-matrices. The uniform distribution is given 
by the invariant measure defined through the relation 

dnp(S) = d^(U'SV), (4) 

which is assumed to be normalized, i.e., 

dfipiS) = 1. (5) 

The matrices U' and V are arbitrary but fixed n x n unitary matrices for (3 = 2, 
V = U' T for (3 = 1, and V = U' R for (3 = A. Equation (0} defines the Circular 
Ensembles, Orthogonal ((3 = 1), Unitary {(3 = 2), and Symplectic {(3 = 4). For these 
ensembles the average of the scattering matrix satisfies (S) = 0. 

When direct processes due to short trajectories exist (prompt responses), the S 
matrix is no longer uniformly distributed. It turns out that in this case the direct 
processes can be characterized by the average (S)(j£ 0), known as the optical S'-matrix. 
In fact, in the maximum-entropy approach developed in the past [Him IB], the probability 
distribution of S is given by 



dP, 



f ) (S)=p$(S)d(l p (S), (6a) 
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where 



7 W (S) = [det(I n -(3) (Syjf ^-m 
\det{I n - S (S) r )\(P n + 2 -V 



is the so-called Poisson kernel. We note that for (S) = the Poisson kernel reduces to 
unity. Therefore, it is clear that the presence of prompt scattering processes makes the 
analysis of the statistical properties of S more difficult. Fortunately, a scattering matrix 
which satisfies the Poisson distribution, Equations (16 ap and (l6H|) . can be transformed 
into a matrix So with uniform distribution by an ad hoc transformation which maps the 
problem with presence of direct processes into one without such processes, namely 

where t c and t' c satisfy [TDl CD] 

t% =I n -(S)HS), (8a) 
t' c t'J = I n -(S)(S)K (8b) 

It has been shown that for one-energy scattering matrices, the Jacobian Jp of the 
transformation given in ([7j) is just the Poisson kernel, (|66|) [31 El El HD] , i-e., 

d^(S ) - J p d^{S) = ldet{In _ s{s y Wn+2 ^ ( 9 ) 

The transformation (J7j) was also used in Ref. [12j to show that the statistical properties 
of the transformed unitary scattering matrices at several energies Sq(Ei), So(E 2 ), ... are 
the ones associated to the problem of scattering in the absence of direct processes. 

We note from (I7j) that for (S) = 0, S reduces to So and a uniform distribution is 
recovered for S, i.e., dP^(S) reduces to the invariant measure dfip^S). Thus, <Q implies 
that if So is uniformly distributed, then the matrix S obtained from (7|) is distributed 
according to the Poisson kernel. 

Although the Poisson kernel distribution has been successfully verified, for example 
in describing several transport properties in quantum dots [31 HI QUI [131 EEH E51 EE] 
and disordered waveguides [17] , there are situations where flux conservation is violated. 
Therefore, the corresponding S matrix becomes non- unitary, i.e., (pD) is not satisfied. 
For instance, power losses (absorption) are unavoidable in experiments on microwave 
cavities and graphs [HI [191 [20] , acoustic resonators [21], and elastic media [22]; in those 
cases the S matrix is sub-unitary. Also, systems with gains (amplification) exist [23J, 
where S is a supra-unitary matrix. 

Recently, several efforts have been made in order to incorporate the information of 
losses or gains (for a review see [21] and [25]; see also [231 [23 [23 12EI (23 EQ])- However, 
few investigations have considered non-unitary scattering matrices S in the presence of 
direct processes. For example, sub-unitary scattering matrices have been considered in 
Refs. [3B E21 E3] • There, the prompt responses come from direct reflections due to an 
imperfect coupling of the antenna to the cavity and the optical scattering matrix (S) 
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was assumed to be a diagonal matrix. In what follows we will consider the general case 
when (S) is a full matrix. 

With the same philosophy of Ref. [8], in this paper we transform a non-unitary 
scattering matrix S, with (S) ^ 0, into a non-unitary scattering matrix So for which 
(So) = 0, by using the transformation (JTj). We show that the Jacobian Jp associated to 
this transformation is given by the square of the Poisson kernel, f|65|) with S replaced 
by S. We consider the cases ofiixn non-unitary scattering matrices in the presence 
(P = 1) and absence (f3 = 2) of time reversal symmetry. Self dual S-matrices ((3 = 4) 
are also studied. 

The paper is organized as follows. The next section is devoted to the invariant 
measure for non- unitary matrices. In Sect. [3] we write So in terms of S using the 
transformation in (jZ|) and the corresponding Jacobian Jp is calculated. The one-channel 
case is presented as a simple example in Sect. HI Finally, we present the conclusions in 
Sect. El 



2. The invariant measure for non-unitary matrices 

Let S be an n x n non-unitary scattering matrix, 

S& + In- (10) 

As for the unitary scattering matrices, we consider the following cases. In the presence of 
time-reversal invariance S = S T and we will refer to this (3 = 1; in the absence of 

time-reversal invariance S has no restriction, and we denote this case by (3 = 2. Finally, 
when S = S R , = 4. 

As in ([3]), any non-unitary matrix S can be parametrized as [9j [33l El] 

S=UpV, (11) 

where p is a real diagonal matrix whose diagonal elements are positive, p ab = p a o~ a b- For 
sub-unitary matrices the eigenvalues of SS^ are real numbers between and 1, while 
they are larger than 1 for supra-unitary matrices. U and V are defined in the same 
way than in (J3j) for each (3 symmetry. Note that for p a = 1, S reduces to a unitary 
S'-matrix and S vanishes for p a = 0. In this sense, the diagonal elements p a determine 
the strength of the absorption or amplification [33] . 

The uniform distribution for an ensemble of S'-matrices is given by the invariant 
measure, defined as in by replacing S by S. In terms of the independent elements 
of S, dfip(S) is given by IB] 

<W(£) = n d(ReS afe )d(ImS\ fe ), for = 1,2. (12) 

{a,b} 

We have used the following notation for the indices for symmetric complex (/3 = 1) and 
complex ((3 = 2) S'-matrices: 

J a < b = 1, . . . ,n for (3 = 1, 
{a ' &}= a,&=l,...,n tor (3 = 2. (13) 
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For self-dual ^-matrices ((3 = 4) with complex quaternion components we have 

d^(S) = J] fl d (iteSff) d (TmS£>) n d d (hnS$) . (14) 

Equations (fT2|) and (fT4|) ensure that (S) = 0, as it should be for a uniform distribution 
of S. 

By differentiation of fTTTj) we obtain 

dS = {Up 1 ' 2 ) 5M {p l ' 2 V) , (15) 

where 

SM = p- 1 ' 2 dU p 1 ' 2 + p- 1 ' 2 dp p- 1 ' 2 + p 1 ' 2 dV p- 1 ' 2 . (16) 

We note that SM is in general a complex matrix ((3 = 2); in addition, it is symmetric 
for (3 = 1 and self-dual for (3 = 4. Alternatively, the invariant measure can be written 

as 

<W(£) = n Re(SM ab )lm(SM ab ), for (3 = 1, 2, (17) 

{a,6} 

and, for /3 = 4, 



<W5) = J] II Re(5MS ) )Im(5MS ) ) nRe(5MS ) )Im(5Mi° 6 ) ), (18) 
where <5M^ is the a-th quaternion component of SM ab (see Appendix A ). 



3. Mapping the scattering problem from presence to absence of direct 
processes 

Let us consider an ensemble of non-uniformly distributed S matrices with average 
(S) 7^ 0. We now apply the transformation given by (J7J) to S matrices in order to 
obtain an ensemble of So matrices with uniform distribution, i.e., 

*° 4 ( s - {S >) l^m* (19) 

where t c and t' c satisfy (similarly to Equations flHaj) and fl85|) ) 

=/„-(£)+<£>, (20a) 
F B e = /„-(5>(5>t. (206) 

Then we look for the Jacobian Jp of this transformation. 
By direct differentiation of S we obtain 

dS = A'dSA, (21) 

with A and A' the complex matrices 

A = (in-iSysY'tl (22a) 



A' = ^(/ n -SW) \ (22&) 
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for (3 = 2. A' = A T for (3 = 1 and A' = A R for (3 = 4. On the other hand, a similar 
expression to ffl5|) is valid for So as well, namely 

dS = (U pl /2 ) 5Mo (pI /2 V ) . (23) 



Substituting Equations (|T5l) and (|23|) into (121]) . we obtain the relation between the 
matrices 8Mq and 5M, 

5M = B'SM B, (24) 

where 5M has the same structure of (fTBl) and 

5 = p^ViVVo 1 ^, (25a) 
B' = Po^U^A'Up 1 / 2 , (25b) 

for p = 2. Again, 5' = S T for /3 = 1 and £' = S fi for /? = 4. We remark that 5M and 
5M are complex matrices, hence B and B' are also complex, in contrast to the case of 
unitary scattering matrices where 5Mq and SM may be taken as real matrices. In the 
following we will restrict ourselves to the cases (3 = 1 and 2, unless explicitly indicated 



otherwise. The case (3 = 4 is considered in Appendix A 



It is convenient to separate the real and imaginary parts of 5M and 5M to obtain 
the Jacobian of the transformation as in (CG 



Re(5M 0ab ) = MB' ac B db )M5M cd ) + ^ Re(iB' ac B db )lm(6M cd ), (26a) 

c,d=l c,d=l 
n n 

lm(5M 0ab ) = MB' ac B db )Re(5M cd ) + £ lm(iB' ac B db )lm(5M cd ). (266) 



Next, we calculate the Jacobian J ab of the transformation which relates the real and 



c,d=l c,d=l 
ab 

imaginary parts of the independent elements of SM with those of 5M as 

Re(6M 0ab )lm(5M 0ab ) = J { ® Re(5M ab )lm(5M ab ). (27) 

To this end, and following Ref. [Sj, we start by assuming that B and B' are diagonal 
matrices, before considering the general case of any complex matrix. 

3.1. A simple example: B and B' diagonal matrices 
Let B and B' 

B a b = K S a b, (28a) 
B' ab = KS ab , (28b) 

where A a 's and A' a 's are complex numbers. Therefore (126 ap and (126 6|) reduce to 

Re(5M 0ab ) = Re(X' a X b ) Re(5M ab ) - lm(X' a X b ) lm(SM ab ), (29a) 
lm(5M 0ab ) = Im(X' a X b ) Re(5M ab ) + Re(X' a X b ) lm(5M ab ). (29b) 

From these two equations, the Jacobian is given by 

Jf = [Re(A' a A fe )] 2 + [Im(A' a A fe )] 2 = |A' a A b | 2 . (30) 
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3.2. The general case 

Any complex matrices B and B' can be written as 

B =QAO, (31a) 

B' = Q'A'0', (31b) 

where Q, Q', O, O' are the most general unitary matrices; we recall that for (3=1, 
B' = B T , hence Q' = T and O' = Q T . Here, A and A' are diagonal matrices whose 



diagonal elements A a and are complex numbers. We substitute (131 ap and (131 fcl) into 
(1241) to obtain 



5M = Q'5M 2 0, (32) 

where we have defined 

5M 2 = A'o"MiA, (33) 

5Mi = 0'5MQ. (34) 

We now calculate the Jacobian of the transformation 5 Mi — > 5M 2 given by (|33|) . 
As before, it is convenient to separate the real and imaginary parts, 

Re(5M 2ab ) = Re(X' a X b ) Re(5M lab ) - lm(X' a X b ) lm(8M Ub ), (35a) 
lm(5M 2ab ) = lm(X' a X b ) Re(5M lab ) + Re(X' a X b ) lm(5M lab ). (35b) 

This is exactly the same transformation as in (129 ap and (129 £>[) for the case of diagonal 
matrices B and B' . Then, the corresponding Jacobian of the transformation in (135 a\i 
and §5M is 

given by (1301) . from which we obtain 
H Re(5M 2ab )lm(5M 2ab ) = J p J] Re(5M lab )lm(5M lab ), (36) 



{a,b} {a,b} 



where 



h = n i 

{a,b} 



(13) 
ab 



{a,b} 



(37) 



In Appendix B we show the explicit calculation of J#. We note, on the other hand, that 
the transformation 5M — > Q'5MO does not change the measure [2]; therefore, from (13"2l 
and (1341) . 5M has the same measure as 5M 2 while 5 Mi has the same measure as 5M. 
Thus, by transitivity, 

II Re(SM Qab )lm(5M 0ab ) = Jp\[ Re(5M ab )lm(5M ab ), (38) 

{a,b} {a,b} 

which means that the Jacobian Jp of the transformation given in (1T9"|) relates the 
invariant measures d/i(So) and dfi(S). Inserting the last result in ([17]) we obtain 

dn(S ) = Jpdfi(S). (39) 

From (1371) for (3=1 and 2, and (1A. 12|) for (3 = 4, Jp can be written in a single 
expression as 

Jp = I (det B') (^+ 2 -/3)/ 2 (det B) 1 2 . (40) 
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Also (ggflD and (1256]) imply that det 5' = det A' and det 5 
expressed as 

J^ = |det {A'A)f n+2 ~ P) . 
Using (122 ap and (122 b\ in the last expression we have 

det t'l det 3 



det A; hence, Jp can be 



(41) 



det (A' A) 



2 ■ 



det [I n -3^S / 
Finally, from (120 ap and (120 6h we can verify that | det t' c \ 



(42) 



det £ c |. Then Jp is given 



by 



det ( I n - (S) (S 



~\ t 



(/3n+2-/3)/2' 



(l(>t ( /„ - S(S^ 



Wn+2-P) 



(43) 



Equation (j4*3l . together with (1391) . is the main result of this work and can be 
interpreted as follows: if a non-unitary scattering matrix So is uniformly distributed 
in the space of non-unitary scattering matrices, then another non-unitary scattering 
matrix S, obtained from Sq through the transformation given by [W\l , is distributed 
according to Jp given by (4^- 111 this sense, Jp is the generalization of the Poisson 



kernel for non-unitary scattering matrices. We show in Appendix C that (1391) together 
with ( 1431) yields the Poisson kernel for unitary matrices. We remark that for unitary 
scattering matrices, the Poisson kernel has been originally obtained in the framework 
of maximum-entropy or Shannon information theory [3]. However, to the best of our 
knowledge, a derivation of the Poisson kernel for non-unitary scattering matrices from 
maximum-entropy arguments is not available. 



4. The one-channel case 



As an example, let us consider a 1 x 1 matrix S which can be parametrized in polar 
form as 



S = VRe 16 , 



(44) 



where R is the reflection coefficient and 6 is the negative of twice the phase shift with 
< R < 1 for sub-unitary matrices and R > 1 for supra-unitary matrices. 

In the sub-unitary uniform distribution for S means that it is distributed 

according to its invariant measure 

d9 



dfipiS) = dR 



2tt' 



A non-uniform distribution of S is constructed from (1451) as 

2tt 



(45) 



(46) 
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For instance, if So is the scattering matrix associated to chaotic cavities with losses in 
the absence of direct processes, p^ (Ro, 6 ) = p ( o\Ro), where p^\R ) is known [21 [26]. 
In the presence of direct processes, according to ( |39|) . with Jg given by ( l4~3i) . we have 



2tt 



S 



\ 



V 



i -s(s x ' 



J 



2tt' 



which, multiplying by po(Ro(R, 9)), can be written as 

/ 



p (Ro(R,9))dR -^ 

Z7T 



1 



s 



2 \ 2 



V 



i-s(s 



J 



Po (R (R,9))dR—. 

ZTY 



Comparing the right hand of side of Equations ( )46l) and ( )48|) . we obtain that 



p w (R,9) 



( 1 



5* 



2 V 



V 



i-sls 



J 



(47) 



(48) 



(49) 



This result, is in agreement with that of Ref. [31], where p^'(R, 9) was verified by 
comparing with experimental measurements in microwave chaotic cavities. 



5. Conclusions 

We have reduced the problem of scattering in the presence of direct processes to the 
case without such processes for n x n non-unitary scattering matrices. We use a 
transformation to map an ensemble of such matrices S with (S) ^ to an ensemble of 
So scattering matrices with (S ) — 0. In our theoretical framework, the direct processes 
are characterized by the average (S). Therefore, S and S describe a system in the 
presence and in the absence of direct processes, respectively. The Jacobian Jp of the 
transformation turns out to be the square of the known Poisson kernel. In our analysis 
we consider general complex, symmetric, and self-dual scattering matrices, in analogy 
to the three basic symmetries in Dyson's scheme /3 = 2, 1, and 4. We have found that 
if So is uniformly distributed in the space of non-unitary scattering matrices, then S 
obtained from So through the transformation ( fl9l) is distributed according to Jp given 
in (l4"3l) . Our study extends known results for the probability distribution of non-unitary 
scattering matrices in the absence of direct process to the case in the presence of such 
processes. As a consequence, from the simplest case of a uniform distribution of S, 
it is possible to obtain more complex distributions, emerging from situations where 
prompt responses are relevant to the scattering problem, as has been illustrated in the 
one-channel case. 
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The invariant measure for (3 = 4 is given by (1181) . For self-dual 5 Mo and 8M matrices, 
(J2H) implies that 

B' = B R . (A.l) 

Using the standard notation for quaternions [2], the elements of SM can be written as 



5M, 



a b 



J25M^e a , a = 0,1, 2, 3, 



a=0 



where SM^ is the projection of SM ab on the quaternion e a , where 
e 



1 




i 




1 




i 


1 


, ei = 


-% 


, e 2 = 


-1 


> e 3 = 


i 



From (1241) . we have 



a,b 



As in Sect. 13.11 we consider first the special case of a diagonal matrix B with 

B 'ab = K$ab- 

Substituting Equations ( 1A.5I) and ( 1A.6I) into ( 1A.41) we obtain 



SMoif = X' i Xj5M-f . 
Hence, the real and imaginary parts are related by 



Re(SM \f 



(a) 



Im(5M 0iJ 



(a)N 



Therefore, we have 

Re(6M %>)lm(5M OiJ 



Re(A^A i )Re(JM i ( ? Q) ) - im^A^im^M^ 
Im(A^A 

i,„fX*4 («) 



Im(XA i )Re(^ a) ) + Re(A^A i )Im(^ c 



A;A i | 2 Re(^ ) )Im(M, 



(A.2) 

(A.3) 
(A.4) 

(A.5) 
(A.6) 

(A.7) 

(A.8) 
(A.9) 

(A.10) 



where |A^Aj| 2 is just the Jacobian of the transformation given by Equations (IA.8I) 
and (1A.9j) . Thus, using ([TBI we obtain 

d/i 4 (^o) = Ja&^a{S) } 



(a.ii; 



where 



J, 



nniwm^ 



(detB 



A2(n+1) 



(detB) 



2(n+l) 



= KdetS^) 2 "" 1 



(detS) 



2n-l I 



(A.12) 



(detB'^detS) 3 
which is the result presented in (j4U|) for (3 = 4. 

As in Sect. 13.21 the general case of a non diagonal matrix B can be reduced to a 
diagonal one by means of Equations (13 lap and (1316j) with Q' = R and O' = Q R . Thus, 
we have the same result for J 4 . 
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Appendix B. Explicit calculation of Eq. ( I37j) 

For (3 — 1, X' a — X a and the products appearing in the right-hand side of Eq. (1371) can 
be written explicitly as 

n n 
a<b a=l b=a 

= (A 1 A 1 )(A 1 A 2 )(A 1 A 3 ) ■•• (AiA n _ 1 )(A 1 A n ) 
x (A 2 A 2 )(A 2 A3) ••• (A 2 A n _i)(A 2 A n ) 

x (A 3 A 3 ) ••• (A 3 A„_i)(A 3 A n ) 

)(A n _iA 

x (A n A n ) 
We can see that Ai appears n + 1 times in the first line only, A 2 once in the first line 
and n times in the second line, A3 appears once in the first and second lines and n — 1 
times in the third line, etc. Then each A a appears n + 1 times in total and the product 
can be written as 

n / n \ n+l 

n AaAf > = n a - +i = n x A = ^ B ^ n+i ^ 

a<b a=l \a=l / 

For (3 = 2, we have 



n /n \ n / n 



n A ^=n ri A ^ =n A "ii A 

a,b=l a=l \b=l J a=l \ b=l 



a=l b=l \a=l / \b=l 

7 ,n . 



= detB' n detB n . (B.3) 
We can summarize the result for any (3 as 

n 

J] A' a A, = det B' {Pn+2 - p)/2 det B^ +2 ~^ 2 , (B.4) 

{a,b} 

which is also valid for (3 = 4 (See Appendix A). 



Appendix C. Reduction to the Poisson kernel for unitary matrices 

We show in this appendix that the original Poisson kernel can be obtained by restricting 
a non-unitary S matrix to be unitary. 

On one hand, we are concerned with unitary scattering matrices So = UqVq and 
5* = UV related by ([7]). Their measures satisfy (Q. Therefore, if the probability density 
distributions of So and S are po(So) and p(S), respectively, po(So) and p(S) are related 
by 

p(S) = po(S ) Jp, or P0 (S ) = ^P-. (C.l) 

Jf3 
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In particular, we assume that p(S) = 5(SS^ — I n ), where S = UpV is a non-unitary 
matrix which is a function of S. Then, So = U p V is a function of So, related to S via 
( TT9l) . This leads to po(So) — S(§o§l — I n ), where So is a function of So- From ( 1C.1I) we 
have that 

5{p -I n ) = %^£. (C.2) 



On the other hand, consider non-unitary scattering matrices, starting with (1391) 
with J/3 given by fH3l) . Here, S is uniformly distributed in the space of non-unitary 
scattering matrices. We restrict to the space of unitary scattering matrices, imposing 
unitarity on Sq by multiplying (139]) by 5(p — I n )- We obtain 

S(po ~ In) dp,{S ) = $(po - h) Jp dfi(S). (C.3) 
Using flCJl, (EJ) yields 

£(po - 4) <Mp ) ^p('S'o) = — j ^ J/3 dfi{p) dp,(S), (C.4) 

where we have written dp(S) = dfi(p) dfi(S). On the left-hand side, the integral with 
respect to p gives one, while the integral with respect to p on the right-hand side 
evaluates to Jp at p — I n . This implies Jg(p = /„) = J| (compare with (H3^) ). 
Finally, we get 

J2 

dfi(So) = -^d/j(S) = Jpdn(S), (C.5) 
where Jg is the Poisson kernel for unitary scattering matrices S. 
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